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ABSTRACT 


In this thesis we have given a dual based procedure for solving the 
transshipment problem . The procedure given by us in this thesis runs in O ( c n^ ) number 
of steps . We found that this procedure reaches 97 % closer to the optimal solution on an 
average . Thus a transshipment problem can reach very close to optimal solution and after 
that it can leave it to conventional 0 ( n^ log (n)) procedures to reach the optimal solution 
. Our contribution is to reach very close to optimal solution in attractive computational 


time 



1 AO 


Chapter No. Description Page No. 

• Ackncwledgement 

• Abstract 

• List of Tables 

• List of Charts 

• Definitiom of Notation 


1. INTRODUCTION 1 

2 LITERATURE REVIEW 5 

3 . PROBLEM FORMULATION & SOLUTION METHOD 

3.1 Problem Formulation 15 

3.2 Solution Method 19 

4. COMPUTATIONAL FORMULATION AND ANALYSIS 

4.1 Introduction 32 

4.2 Small Size Transshipment Problems 35 

4.3 Medium Size Transshipment Problems 40 

4.4 Large Size Transshipment Problems 42 

4 5 Remark on no. of Iterations taken by our Heuristic 49 

4.6 Running the Heuristic in 0( c n^ ) times. 50 

5. CONCLUSIONS AND DIRECTIONS FOR FUTURE 

RESEARCH 57 

REFERENCES 58 

APPENDICES . 

Appendix I’ Solved Numerical Examples. 61 

Appendix II : Pascal Program Listing . 71 

CHARTS : 


• Charts showing percentage closeness to the optimal result. 

• Charts showing comparision between no. of iterations taken by the heuristic 
procedure and CPLEX. 



LIST OF TABLES 


Problem Size 

Description 

Number of 

Problems 

Tested 

Page 

No 

3x3x3 

Small Size Problems 

15 

36 

4x4x4 

Small Size Problems 

15 

37 

5x5x5 

Small Size Problems 

10 

38 

7x7x7 

Small Size Problems 

10 

39 

10x10x10 

to 

75 X 75 X 75 

Medium Size Problems 

10 

41 

100 X 100 X 

100 

Large Size Problems 

50 

43 

150 X 150 X 

150 

Large Size Problems 

50 

45 

200 X 200 X 

200 

Large Size Problems 

50 

47 

100 X 100 X 

100 

Results from 0(c ) 

Solution Procedure 

50 

51 

150 X 150 X 

150 

Results from 0(c ) 

Solution Procedure 

50 

53 

200 X 200 X 

200 

Results from 0(c ) 

Solution Procedure 

50 

55 























































LIST OF CHARTS 


Charts Showing Percentage Closeness 


Problem Size 

Page No. 

3x3x3 

83 

4x4x4 

83 

5x5x5 

84 

7x7x7 

84 

100 X 100 X 100 

85 

150x 150x 150 

86 

200 X 200 X 200 

87 


Charts Showing Comaprisions of Iterations Taken 
by Heuristic and CPLEX 


Problem Size 

Page No. 

3x3x3 

88 

4x4x4 

88 

5x5x5 

89 

7x7x7 

89 

100 X 100 X 100 

90 

150 X 150 X 150 

91 

200 X 200 X 200 

92 





DEFINITIONS OF NOTATION 


Notation 


Definition 



XI, j 



C2a- 



Ik 


WJ 


C22k 



Number of Plants 


Number of Intermediate PointsAVarehouses 


Number of Markets 


The Quantity Recieved as a Fraction of Total Market 
Demand from Plant i to Intermediate Point j 


The Quantity Recieved as a Fraction of Total Market 
Demand from Intermediate Point j to Market k 


The Absolute Value of Demand at Market k 


The Demand at Market k as a Fraction of Total Market 
Demand 


K 

Cost of Shipping 2] Dk units from Plant i to 

Ar=l 

Intermediate Point j 


K 

Cost of Shipping ^ Dk units from Intermediate Point j 

k=\ 

to Market k 


The Absolute Value of the Plant Capacity at Plant 
Location i 


The Capacity of the Plant i as a Fraction of total Market 
Demand 


min (Clij ) Vi = 1..I 

J 


min(C2jk) Vk = l..K 

J 


: C2k — C2i k 


Set of Wj’s that can be Increased 


min C2jk Vk 


The Dual Variables Associated with Constraints 
(2),(3),(4),(5),(6) respectively . 


The Sum of Supply, Intermediate and Demand Points 



































CHAPTER 1 


INTRODUCTION 

The uncapacitated transshipment problem is a special case of minimum cost flow 
problem which includes all the general features of minimum cost flow problem except for 
not having (finite) arc capacities . In Transshipment problem , flow takes place in stages, 
fi'om plant goods move to warehouses and fi-om warehouses they move to the markets. 

The Transshipment problem is an example of linear network optimization 
and is now a standard application for industrial firms having several manufacturing plants , 
warehouses , sales territories and distribution outlets . In this instance , the strategic 
decisions involve selecting transportation routes so as to allocate the production of 
various plants to several warehouses, or terminal points. 

A large class of network flow problems centre around the shipping and 
distribution problems . This problem is best described in terms of shipments from plants to 
warehouses and fi'om warehouses to the markets . 

Network flow is an important area in linear programming and is known for 
its novel theoretical properties . Most network flow problems are linear programming 
problems and their special structure often allow us to develop highly eflBcient algorithms 
that do not seem possible for the general linear programming problems . This eflSciency in 
solving network flow problems has motivated the practitioners to identify the class of 
problems that can be modelled as network flow problems . Many successes have been 
made at this front and a wide variety of practical problems not only in transportation and 
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communication , but also in facility location , production planning , project management , 
machine loading , operation scheduling and cash management can be formulated and 
solved as network flow problems . Recent developments in solution methodology , 
implementation technology and improved computer software have enabled the network 
flow problems to be handled much more effectively and conveniently . 

Network flow problems are found in various types . Some well known types of network 
flow problem are : 

1 The Shortest Path Problem 
2. The Maximum Flow Problem 
3 The Transportation Problem 

4. The Transshipment Problem 

5. The Assignment Problem 

6. The Minimum Cost Flow Problem 

In the above mentioned types of network flow problems , the first five 
types (i.e 1-5 ) are special cases of the minimum cost flow problems and can be solved 
even more efficiently than the minimum cost flow problems . The minimum cost flow 
problems is however the most popular network flow problem . Some of the reason for its 
popularity include ' 

• Simplicity and intuitive appeal of the problem as well as its solution procedures ; 

• Integer valued optimum solution ; 

• Little computer memory requirement ; 

• EfiBciency of computer codes in solving even the large sized problems ; 

• The ability to formulate several complex problems as a sequence of the minimum cost 


flow problems . 


Due to these reasons , the minimum cost flow problem has been 
extensively researched in the literature . Researchers have developed numerous algorithms 
. In the past , polynomial algorithms for the minimum cost flow problems have been 
vigorously pursued and several new algorithms have been developed . Among these 
contributions are the cost scaling algorithm due to Goldberg and Taijan [1987] and the 
double scaling algorithm due to Ahuja, Goldberg, Orlin and Taijan [1988] . These 
algorithms have very attractive worst case time bounds but their empirical performance is 
unclear . 

In most of the above mentioned cases , the well known primal approaches 
solve the uncapacitated transshipment problem in O (n^ log(n)) times . In this paper we 
tried to give a dual based procedure which takes O (c n ^ ) times to solve the 
transshipment problem , where n is the number of nodes in the network . Though the 
solution given by us in this paper does not give an optimal solution, but covers most of its 
path towards optimality in O ( c n^ ) time and thereafter the conventional dual based 
algorithms can take over to reach the optimal solution 
The organization of the thesis is as follows : 

In chapter two , we have given relevant literature review which highlight the 
various literature available in the field of the minimum cost flow problems inclusive of the 
best available algorithm at present . The problem formulation and the heuristic solution 
procedure to solve the transshipment problem is given in chapter three . In the following 
chapter four, we have attempted to provide our computational experience and analysis 
obtained during the testing of the heuristic on small as well as large sized randomly 
generated problems . In chapter five we wind up by giving conclusion and directions for 
future research in this field . Appendix I contains few illustrated numerical examples to 
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provide illustration of the solution procedure given by us in this thesis . Appendix II 
contains the computer program developed by us to carry out the step by step heuristic 
solution procedure . We also attached various graphs to illustrate our computationally 
efficient 0( cn^ ) procedure 

Our solutions are within 98 % of the optimal solution . Our solution can be 
used as an advanced start by the conventional dual based approaches for the 
transshipment problem One more interesting feature of our thesis is that an 0( c n^ ) 
procedure developed by us took us to 96.5 % closer to the optimal solution . 



CHAPTER 2 


LITERATURE REVIEW 

The minimum cost flow problem is an important network flow problem . 
The transshipment problem , a special case of minimum cost flow problem , was posed 
and solved by A. Orden [1956] , Koopmans [1947] and Kantorovnch [1939] . 

Minty [1960] and Fulkerson [1961] independently discovered the out-of- 
kilter algorithm . Aashtiani and Magnanti [1976] have described an efficient 
implementation of this algorithm . The cycle-canceling algorithm is credited to 
Klien[1967] Three special implementations of the cycle-canceling algorithm run in 
polynomial time ; the first , due to Barahona and Tardos [1989] ( which , in turn , 
modifies an algorithm by Weintraub [1974] ) , augments flow along ( negative ) cycles 
with the maximum possible improvements ; the second , due to Goldberg and 
Taijan[1988] , augments flow along minimum mean cost (negative) ; the third , due to 
Wallachar and Zimmerman [1991] , augments flow along minimum ratio cycles . 

Helgason and Kennington [1977] and Armstrong , Klingman and 
Whitman[1980] describe the specialization of the linear programming dual simplex 
algorithm for the minimum cost flow problem . Each of these algorithms perform 
iterations that can (apparently) not be polynomially bounded. Zadeh [1973 a] describes 
one such example on which each of several algorithms - the primal simplex algorithm with 
Dantzig’s pivot rule , the dual simplex algorithm , the negative cycle algorithm , the 
successive shortest path algorithm , the primal dual algorithm , and the out-of-kilter 
algorithm - performs an exponential number of iterations . Zadeh [1973 b] has also 
described more pathological examples for network algorithms . 



The network simplex algorithm and its practical implementations have 
been most popular with operations researchers . Johnson [1966] suggested the first tree 
manipulating data structure for implementing the simplex algorithm . The first 
implementations using the ideas , due to Srinivasan and Thompson [1973] , and Glover , 
Klingman and Napier [1974], significantly reduced the running time of the simplex 
algorithm Glover, Klingman and Stutz [1974], Bradley, Brown and Graves [1977] and 
Barr, Glover and Klingman [1979] subsequently discovered improved data structures . 
The book of Kennington and Helgason [1980] is an excellent source for references and 
background material concerning these developments . 

The relaxation algorithms proposed by Bertsekas and his associates are 
other attractive algorithms for solving the minimum cost flow problems and its 
generalization . For the minimum cost flow problem , this algorithm maintains a pseudo 
flow satisfying the optimality conditions The algorithm proceeds by either 

• augmenting flow fi'om an excess node to a deficit node along a path consisting of arcs 
with reduced cost, or 

• changing the potentials of subset of nodes . 

In the later case , it resets flows on some arcs to their lower or upper 
bounds so as to satisfy the optimality conditions ; however , this flow assignment might 
change the excesses and deficits at nodes . The algorithm operates so that each change in 
the node potentials increases the dual objective fimction value and when it finally 
determines the optimal dual objective function value , it has also obtained the optimal 
primal solution . This relaxation algorithm has exhibited nice empirical behaviour . 
Bertsekas [1985] suggested the relaxation algorithm for the minimum cost flow problem 
(with integer data ) . Bertsekas and Tseng [1985] extended this approach for the 


minimum cost flow problem with real data and the generalized minimum cost flow 
problem . 

Ahuja , Batra , Gupta and Punnen [1996] address the problem of allocating 
a given budget to increase the capacities of arcs in a transshipment network to minimum 
cost flow in the network, the capacity expansion costs of arcs are assumed to be piece 
wise linear convex functions . They used the properties of the optimum solution to 
convert this problem into a parametric network flow problem . The concept of optimum 
basis structure is used which allows to consider piece wise linear convex functions 
without introducing additional arcs . The resulting algorithm yields an optimum solution 
of the capacity expansion problem for all budget levels less than or equal to the given 
budget . For integer data this algorithm performs almost all computations in integer. 

A number of empirical studies have extensively tested minimum cost flow 
algorithms for wide variety of network structures , data distributions , and problem size. 
The most common problem generator is NETGEN , due to Klingman , Napier and 
Stutz[1974] , which is capable of generating assignment, capacitated or uncapacitated 
transportation / transshipment and minimum cost flow problems Glover , Kamey and 
Klingman [1974] and Aashtiani and Magnanti [1976] have tested the primal dual and out- 
of-kilter algorithms . Helgason and Kennington [1977] and Armstrong , Klingman and 
Whitman [1980] have reported on extensive studies of the dual simplex algorithm . The 
primal simplex algorithm has been a subject of more rigorous investigation ; studies 
conducted by Glover , Kamey , Klingman and Napier [1977] , Glover, Kamey and 
Klingman [1974] , Bradley , Brown and Graves [1977] , Mulvey [1978 b] , Grigoriadis 
and Hsu [1979] and Grigoriadis [1986] are noteworthy . Bertsekas and Tseng [1988] 
have presented computational results for the relaxation algorithm . 



Bertsekas and Tseng [1988 b] developed the relaxation algorithm and 
conducted extensive computational investigations of it . A FORTRAN code for the 
relaxation algorithm appears in Bertsekas and Tseng [1988 a] . Their study and those 
conducted by Grigoriadis [1986] and Kennington and Wang [1990] indicate that the 
relaxation algorithm and the network simplex algorithm are the two fastest available 
algorithms for solvdng the minimum cost flow problem in practice . When the 
supplies/demands at nodes are relatively small , the successive shortest path algorithm is 
the fastest algorithm . Previous computational studies conducted by Glover , Kamey and 
Klingman [1974] and Bradley , Brown and Graves [1977] have indicate that the network 
simplex algorithm is consistently superior to the primal-dual and out-of-kilter algorithms. 

Polynomial Time Algorithm 

In the recent past , researchers have actively pursued the design of fast 
(weakly) polynomial and strong polynomial time algorithms for the minimum cost flow 
problems . Recall that a algorithm is strongly polynomial-time if its running time is 
polynomial in the number of nodes and arcs , and does not involve terms containing 
logarithms of C or U . The table given below summarizes these theoretical developments 
in solving the minimum cost flow problem . The table also reports running time for 
networks with ‘n’ nodes and ‘m’ arcs , m' of which are capacitated . It assumes the 
integral cost coefficients are bounded in absolute value by U . the terms S(.) is the running 
time for the shortest path problem and the term M(.) represents the corresponding 
running time to solve a maximum flow problem. 



weak Polynomial- 1 ime Aigontnms 


S.No. 

Discoverers 

Year 

Running time 

1 

Edmonds and Karp 

1972 

0((n+m) log U S(n , m, nC)) 

2 

Rock 

1980 

0((n+m) log U S(n , m, nC)) 

3 

Rock 

1980 

0((nlogCM(n,m,U)) 

4 

Bland and Jensen 

1985 

0((m log C M(n , m , U)) 

5 

Goldberg and Taijan 

1987 

0((nm log (n^/m) log (nC)) 

6 

Goldberg and Tarjan 

1988 

0((nm log (n) log (nC)) 

7 

Ahuja , Goldberg, Orlin and 

Taijan 

1988 

0((nm log( log U) log (nC)) 

































Strongly Polynomial-Time Algorithms 


S. No. 

Discoverers 

Year 

Running Time 

1 

Tardos 

1985 

0(m^) 

2 

Orlin 

1984 

0((n+my log n S(n,m)) 

3 

Fujishige 

1986 

0((n+my log n S(n,m)) 

4 

Gain and Tardos 

1986 

0((n^ log n S(n,m)) 

5 

Goldberg and Taijan 

1987 

0((nm^ log n log (nVm)) 

6 

Goldberg and Tarjan 

1988 

0(nm^ log^ n) 

7 

Orlin 

1988 

0((n+m') log n S(n,m)) 


For the sake of comparing the polynomial and strongly polynomial-time 
algorithms, we make the assumption the C= 0(n'‘) and U=0(n‘‘) for the same constant ‘k’ 
. Under this assumption the best bounds for the shortest path and maximum flow problem 
are : 


Polynomial Time Bounds 

Discoverers 

S(n,m,c)=0(Min(m+nVlogC),(m log log C)) 

Ahuja,Mehlhom,Orlin and 


Taijan [1990] 

S(n,m,c)=0(Min(m+n-y/logC),(m log log C)) 

Van Emde Baos, Kaas and 


Zijistra [1977] 


Ahuja,Orlin, andTauan[1989] 


m /V 











































Strongly polynomial-Time Bounds : 


S(n,m) = 0(m+n log n) 

Fredman and Tarjan [1984] 

M(n,m)=0(miii(iim+n^'^ ^ ,mn log n) 

King, Rao, Taijan [1991] 


where e is any fixed constant. 

Now we give some well known scaling algorithms used in the study of 
minimum cost flow problems along vdth their running times ; 


Algorithm 

Running Time 

Capacity Scaling Algorithm 

0((m log U)(m+n log n)) 

Cost Scaling Algorithm 

0(n^ log (nC)) 

Double Scaling Algorithm 

0(mn log U log (nC)) 

Minimum Mean Cycle-Canceling Algorithm 

0(nW log n) 

Repeated Capacity Scaling Algorithm 

0((m^ log n)(m+n log n)) 

Enhanced Capacity Scaling Algorithm 

0((m log n)(m+ n log n)) 


Using capacity and right-hand-side scaling, Edmonds and Karp[1972] 
developed the first (weakly) polynomial - time algorithm for the minimum cost flow 
problem , known as the RHS-scaling algorithm for the minimum cost flow problem . A 
variant of the Edmonds-Karp algorithm , was suggested by Orlin [1988] . the scaling 
technique did not initially capture the interest of many researchers , since they regarded as 
having little practical utility. However , researchers gradually recognized that the scaling 
technique has great theoretical value as well as potential practical significance 

























Rock[1980] developed two different bit-scaling algorithms for the minimum cost flow 
problems , one using scaling and the other using cost scaling . This cost scaling algorithm 
reduces the minimum cost flow problem to a sequence of 0(n log C) maximum flow 
problems Bland and Jensen [1985] independently discovered a similar cost scaling 
algorithm 

The pseudo flow push algorithms for the minimum cost flow problem use 
the concept of approximate optimality, introduced independently by Bertsekas [1979] and 
Tardos [1985] . Bertsekas [1986] developed the first pseudo flow push algorithm This 
algorithm was pseudopolynomial-time . Goldberg and Tarjan [1979] used a scaling 
technique on a variant of this algorithm to obtain a generic pseudo flow push algorithm 
Taijan[1984] proposed a wave algorithm for the maximum flow problem, the wave 
algorithm for the minimum cost flow problem was developed independently by Goldberg 
and Taijan [1987] and Bertsekas and Eckstein[1988] , relies upon the similar ideas. Using 
a dynamic tree data structure in generic pseudo flow push algorithm, Goldberg and Taijan 
[1987] obtained a computational time bound of 0(nm log n log(nC)) . They also showed 
that the minimum cost flow problem can be solved using 0(n log (nC)) blocking flow 
computations. Using both finger tree and dynamic tree data structures, Goldberg and 
Taijan [1982] obtained a 0(nm log(n^/m) log(nC)) bound for the wave algorithm, this 
algorithm was subsequently improved to 0(nm log U log (nC)) by Ahuja, Goldberg, Orlin 
and Taijan[1982] . 

The first strongly polynomial-time minimum cost flow problem is due to 
Tardos[1985]. Several researchers including Orlin[1984], Fujishige[1986], Gail and 
Tardos[1986] and Orlin[1988] provided subsequent improvements in the running time 
Goldberg and Taijan [1988 a] obtained another strongly polynomial-time algorithm by 



slightly modifying their pseudo flow push algorithm. Goldberg and Taijan [1988 b] also 
show that their algorithm that proceeds by canceling minimum mean cycles is strongly 
polynomial time Currently, the fastest strongly polynomial-time algorithm is due to 
Orlin[1988] 

Some additional polynomial-time minimum cost flow algorithms include 
(1) a triple scaling algorithm due to Gabow and Taijan [1989 a], (2) a special 
implementation of the cycle canceling algorithm developed by Garahona and 
Tardos[1989], and (3) (its dual approach) a cut canceling algorithm proposed by Ervolina 
and McCormick [1990 a] . 

Interior point linear programming algorithms are another source of 
polynomial-time algorithms for the minimum cost flow problem. Among these, the fastest 
available algorithm, due to Vaidya[1989] , solves the minimum cost flow problem in 
0(n^ ^ ^mK ) time, with K= log n + log C + log U. 

Currently the best available time bound for the minimum cost flow problem 
is 0(min { nm log (n2/m) log (nC)) , nm ( log log U) log (nC), (m log n)(m+n log n)}) ; 
the three bounds in this expression are, respectively, due to Goldberg and Taijan[1987], 
Ahuja, Goldberg, Orlin and Taijan [1992] and Orlin[1988] . At this time the research 
community has yet to develop sufficient evidence to fully assess the computational worth 
of scaling and interior point linear programming algorithms for the minimum cost flow 
problem 

Ali, Padman and Thiagarajan [1989] reported that the number of pivots in 
implementations of the dual algorithm is smaller that in the primal simplex algorithm. This 
feature has obscured the true performance characteristics of the dual algorithms Sharma 



and Sharma [1998] also gave a dual based procedure to solve the simple transportation 
problem in O(cn^) time . 

In the present work we also tried to extend the dual feature as given by 
Sharma and Sharma[1998] to the transshipment problem. We , in this thesis , have 
developed a heuristic solution procedure for solving the transshipment problem in 0(c n^ ) 
time, the details of which is given in Chapter three . 



CHAPTER 3 


PROBLEM FORMULATION 
AND 

SOLUTION METHOD 


In this chapter we give the problem to be solved in this paper and the dual 
based solution procedure to solve the problem We start by defining the constants used 
by us in the problem and subsequently we have given the complete problem formulation 
and the dual based solution procedure 

3.1 PROBLEM FORMULATION 


3.1.1 Constants of the Problem : 

In this paper we used the index i for supply nodes , index k for demand 
nodes and index j for the intermediate / transshipment nodes . 

D k is the demand at demand node k and d t is the demand at demand node k as a 
fraction of total demand Hence, 

D* 

dk =■ 

K 

k=\ 

and obviously we have, 

K 

k=\ 

B, is the supply available at source node i and b j is the supply available as a 


fraction of total market demand i.e. , 



and we assume that , 


t*. = 1 

1 = 1 

Cl ,j is the transportation cost from source node i to intermediate node j for 

K 

transporting the entire market demand ^ Dk 

*=i 

C2 j k is the transportation cost from intermediate node j to demand node k for 

K 

transporting the entire market demand ^ Dk 

*=i 


[In this paper we assume that total supply from all the sources is equal to the total 
demand of all the demand nodes . ] 

3.1.2 Decision Variables : 

XI ,j is the quantity received as a fraction of total demand at intermediate node j 
from source node i 

X2jk is the quantity received as a fraction of total demand at demand node k from 


intermediate node j. 


3.1.3 Problem Formulation : 


We state the problem to be solved in this paper as foUows ■ 



Problem P : 


Minimize Z 

Oj) 

Cl.j . XI, J + 

2 C2jk.X2jk 


Subject to • 

z 

Oj) 

XI, J =1 ViJ 


(2) 


j 

-z 

;=i 

X1,J > -b. 

V i 

(3) 


z 

X2jk = 1 

Vj,k 

(4) 


j 

-z 

J=1 

X2jt > -dk 

V k 

(5) 


1 

z 

Xl,j - 2 X2jk 

= 0 

(6) 


1=1 t=i 


0 ) 







Xlij >0 ViJ (7) 

X2jk > 0 Vj,k (8) 


In the model defined above , we assumed that the total supply at all supply nodes is equal 
to the total demand at all the demand nodes ,which makes it a valid formulation of the 
transshipment problem 

3.1.4 Dual of the Problem P : 

Now we define the dual of the problem P defined above . In the dual of the 
problem , we associate the dual variables Vlo , z, , V2o , Vk , and Wj with constraints (2) , 
(3) , (4) , (5) and (6) respectively . 

The dual of the problem model P (1) - (8) is as follows . 


Problem DP : 



1 

Maximize Zjuai = (Vlo- ^ z, .b, ) + 

j=l 

(V2o-|; dk.Vk) 

( 9 ) 

Subject to : 



Vlo-z> < Clij - Wj 

V ij 

(10) 

V2o -Vk < C2jk + Wj 

V J , k 

(11) 

Zi > 0 Vi 


(12) 

Vk ^ ^0 V k 


(13) 

Vlo,V2o,Wj Vj 

unrestricted. 

(14) 



3.2 SOLUTION PROCEDURE 


The solution procedure is based on the dual of the problem P i.e. the problem DP 
and it aims to increase the value of Z dual as much as possible . We start by defining a 
special case of the problem i e the case when all Wj’s are zero’s 
3.2.1 Special Case of the Problem 

In this special case of the problem DP i.e. when all Wj’s are zero the 
problem DP has the following form ■ 

Problem DPI ; 


J 

Maximize (Vlo - ^ z, . 

1=1 

b.) + (V2o-X dk.Vk) 

*=1 

(15) 

Subject to 



Vlo-z, < Cl.j 

V i,j 

(16) 

V2o - V k ^ C2jk 

V j , k 

(17) 

z, > 0 

V i 

(18) 

Vk ^ 0 

V k 

(19) 

Vlo , V2o 

V j unrestricted. 

(20) 


Problem DPI separates into two independent problems which are identical 
to problem DP2 given below - 

Problem DP2 

K 

Maximize Vo dk Vk (21) 

i=l 

Subject to, 

Vo - Vk < Cik V i,k 

Vk > 0 


Vo unrestricted 


( 22 ) 

(23) 

(24) 



The two independent problems are DP 11 and DP 12 which are given below 


Problem DPI 1 . 


j 

Maximize (Vlo - ^ z> • b, ) 

1=1 


Subject to, 


Problem DP 12 


( 1 6), ( 1 8) and Vlo unrestricted 


K 

Maximize (V2o dt . vt ) 

it=i 

Subject to, 

(17) , (19) and V2o unrestricted 


Sharma[26] has given an efBcient algorithm to optimally solve the problem 
DP2 and it is reproduced here for the sake of completeness 
3.2.1. 1 Algorithm for solving the Problem DP2 : 

Step 1 Find Ck* = min(C, k) for all k = 1 . K and remove all the redundant 

1 

constraints (of the equation (22)) in problem DP2 
In case of a tie , only one equation is retained and all the others are 
eliminated. Then the problem DP2 becomes the following 
Problem DP21 

K 

Maximize (Vo-^ dk-Vk) (25) 

t=i 


V k 


Subject to • 

Vo -Vk Ck* 


(26) 


Vk ^ 0 V k 

(27) 

Vo unrestricted 

(28) 


Step 2 : We sort the values of C k* in an increasing order and reindex such that 

Q’< C^i*forallr=l..K-l 

K 

Step 3 ; Since ^ ct = 1 

*=i 

we set, Vk = Vo - Ck* for all k. 

It has been shown in Sharma[1998] and Shamia[1996] that above 


algorithm solves the problem DP21 optimally 



3.2.2 Definitions : 


Now , in this section , we give some of the definitions which are frequently used in 
the solution procedure . 


Cl.* ; 

min (Cl ,j ) 

j 

V i,j 

m, 

U ci>ci„} 

ViJ 

C2k* • 

min(C2,i) 

J 

V j , k 

It 

{j - C2k*= C2,k} 

V j,k 


WJ : Set of Wj ’s to increase the value of the objective value Z dual 


min C2jk V k 

j £{1 J}-WJ 


C22k 


3.2.3 Objective : 

We want to increase a collection of Wj s . j e WJ so that objective value (Zduai ) 
of the problem DP [(9)-(14)] increases 

If we increase Wj’s then Cl,* for the problem DP 11 decrease and objective 
function value of problem DP 11 decreases but at the same time C2k* for the problem 
DP12 increases and hence objective function value of problem DP12 increases. We have 
to increase Wj’s ■ j gWJ in such a manner as to strike the balance These are the ideas of 
all the heuristics presented below 

3.2.4 Heuristic H 1 : 

In this section we give the heuristic which calculates the net benefit obtained by 
unit increase of the variables W/s . j e WJ . The execution of this heuristic tells whether 
it will beneficial to increase the combination of Wj’s . j e WJ or not .Only if the net 
benefit (as calculated using this heuristic) is greater than zero, we increase the currently 
obtained combination of Wj’s . j e WJ otherwise we go to find the next combination of 
Wj’s ; j € WJ to be increased The Heuristic is defined as follows : 

Given a set of WJ , we determine if it is beneficial to increase Wj’ s : j e WJ 
so that the objective value of the problem DP [(9)-(14)] increases 
Step 1 : Set benefit =0. 

Compute Cl,* , m, Vi [as defined above] 

Compute C2k* ,1k V k [as defined above] 


Step 2 : 


for k = 1 to K do 

if Ik c WJthen 

Set benefit = benefit + dk 


Step 3 Set loss = 0 

Step 4 for j = 1 to J 

if ( j e m , ) and ( j e WJ) then 
set loss = loss + b, 

[ Note ■ b ,’s to considered only once for each i . ] 

Step 5 Compute Net_Benefit = benefit - loss 

If ( Net_Benefit > 0) then 

It is possible to improve the objective functional value of the problem DP 
by increasing the value of W, s . j e WJ by at least one unit (We assume that all costs are 
integer. ) 

Result 1 : The complexity of the heuristic HI in 0(n) times 

Proof : In the heuristic HI defined above the steps 1 is executed in 0(I+K) steps 

while step 2 in done in 0(K) steps and step 4 is done in 0(J) steps . Thus the Heuristic HI 
runs in 0(n) number of steps where n = I + J + K ,ie the sum of supply, demand and 
intermediate points 

Now we give two heuristics below to determine the set of WJ (as defined 


above). 


3.2.5 Heuristic H2 : 


Using this heuristic , we find the combination of W/ s • j e WJ which is 
required to increase the value so that the value of the objective function value Z^uai can be 
improved In this heuristic , we calculate the combination of Wj’ s . j e WJ sequentially 
i.e. taking each Ik at a time as the set WJ . The detailed heuristic is given below in easy 
steps 

Step 1 ; Compute Cl,*,m, Vi [as defined above] 

Compute C2k* ,1k, C22 k V k [as defined above] 

Step 2 • for k = 1 to K do 
Begin 

SetWJ = lk 

Compute Net_Benefit by using Heuristic HI 
If (NetJBenefit > 0) then 

Stop. We have found a collection of Wj s j g WJ 
which can be increase in value to improve the 
objective function value of the problem DP. 

Stop. 

End. 

Result 2 : The Heuristic H2 runs in O(n^) number of steps 

Proof : In the heuristic H2 abpve , step 1 is executed at the most 0(I+K ) steps 

and step 2 is executed at the most 0(K^) times. Hence heuristic H2 is 
executed for at the most O(n^) number of steps where n = I+J+K 


Now we define another heuristic for finding the set WJ 
3.2.6 Heuristic H3 : 

Now in this section we give one more heuristic which can be used to find 
the combination of W/ s : j € WJ to improve the objective function value of Zdmi ■ In 
this heuristic , we find the combination of Wj’ s : j e WJ in a different way. Here we take 
one Ik and find out the common elements of this particular Ik with all others Ik’s . After this 
we include all such common element in the set WJ . For a better understanding of this 
heuristic we give below the complete heuristic in steps below ; 

Step! • Compute CI.*,m, Vi [as defined above] 

Compute C2k’ ,1k, C22 k V k [as defined above] 

Step 2 . For kl = 1 to K do 
Begin 

WJ = lki 

Fork2 = (kl+l)toKdo 

If 1 kl and 1 k 2 both have some common elements then 
Set WJ = WJ + 1 k 2 

Compute Net_Benefit by using Heuristic HI 
If (Net_Benefit > 0) then 

Stop ' We have found a collection of Wj s ; j e WJ 
which can be increase in value to improve the 
objective function value of the problem DP. 


Stop 


End 


Result 3 : The heuristic H3 runs in O(n^) number of steps . 

Proof : In the above stated Heuristic H3, step 1 is executed for at the most 0(I+K) 

number of steps whereas the step 2 executes for at the most 
0(Kl*K2*0(n)) number of steps Hence the complete heuristic H3 runs in 
O(n^) times where n =I+J+K 

After defining the heuristics for determining the set WJ and calculating the net 
benefit which enable us as to which particular combination of Wj’ s j e WJ are to be 
increased finally we require upto how many units we can really increase the values of the 
Wj’ s • j € WJ to improve the objective function value of Zduai • This task is given by us 
in heuristic H4 defined in the coming section 

3.2.7 Heuristic H4 : 

This heuristic calculates the extent of increase of the combination of Wj’ s 
j e WJ to improve the value of objective function Zduai • To calculate this, we first find 
out that upto which extent we can increase the values of W j’ s : j e WJ . For this 
purpose we find extent of increase for each k as the difference between the smallest and 
the second smallest value of C2jk of that particular k . After this, we take the minimum of 
all such values obtained We also calculate that upto which extent we can decrease the 
values of Wj’ s : j e WJ as it also effects the improvement of the solution To find the 
extent of decrease , we first find out the extent of decrease for each value of matrix Cl 
such that j G WJ . Then we find for each i , the extent of decrease possible in a specified 


manner [given below] and take the minimum of all such values to finally obtain the value 
upto which we can decrease the values of W j’ s : j e WJ . Now after calculating the 
value of possible increase and possible decrease , we finally take the minimum of the two 
which gives us finally the value with which we can update the cost coefficients This 
whole procedure is described below in detailed steps 

Step 1 : Given a set of WJ. 

Set Extent_of_Increase k = oo V k 

Step 2 ; Set Extent_of_Increase k = C22 k - C2k* 

Step 3 ; Extent of_Increase2 = min [ Extent_of_Increase k ] 

k 

Step 4 : Extent_of_Decrease ij { j € WJ } 

= Cl.j-Cl.* if Cl ij -Cli*>0 

= 00 if Cl ,j - Cl i* = 0 

Extent_of_Decrease , 

= min {Extent of Decrease ,j ) VExtent_ofJDecrease ,j qt oo 

j ewj 

= 00 If 3 Extent_Decrease ,j = oo 

Step 5 : Extent_of_Decreasel = min {Extent_of_Decrease i) 

1 

Step 6 : Extent_of_lncrease = min{ Extent_of_Increase2, Extent_ofJDecreasel } 


Stop 


Result 4 : 


The Heuristic H4 runs in 0(n^) times . 

Proof : In the above defined heuristic, Steps 1,2,3 runs in 0(K) number of steps 

Step 5 runs in 0(1) number of steps .Step 4 is executed for 0( I*J ) 
times. Hence the computational complexity of the heuristic H4 is O(n^) 
times where n = I+J+K . 

Finally we give a complete Heuristic H5 which we used to increase the 
objective function value of the problem DP using the heuristics HI - H4. 

3.2.8 Heuristic H5 : 

In this section we describe the complete solution procedure which uses all 
the above defined heuristic to obtain the best solution . This whole procedure is given in 
the following lines in very easy steps 

Step 1 ■ Compute Cl,*,m, Vi [as defined above] 

Compute C2k* ,1k V k [as defined above] 

Step 2 Find set WJ using heuristic H2 

Step 3 : Compute Net_Benefit for the set WJ using HI . 

If (Net_Benefit > 0) then 
Go to Step 6 ' 

Otherwise Go to Step 4. 

Step 4 ; Find set WJ using Heuristic H3. 

Steps : Compute Net_Benefit for the set WJ using HI 


Step 6 . 
step 7 . 

Step 8 ; 
Step 9 ; 

Result 5 : 
Proof : 


If (NetJBenefit > 0) then 
Go to Step 6 
Otherwise Go to Step 9 

Calculate Extent_of_Increase using heuristic H4. 

Update costs coefficients Cl ij and C2jk as follows : 

Cl ,j = Cl ,j - Extent_Of_Increase. V j e WJ 

C2j k = C2jk + Extent_Of_lncrease. V j g WJ 

Go to Step 1 

Stop. Print the best solution obtained so far. 

Thus with the heuristic H5 , we found the best solution 
Heuristic H5 runs in 0( c n^) number of steps . 

In heuristic H5 defined above, step (1) runs in 0( I+K) number of steps 
Step (2) [H2] runs in 0(n^ ) times Step (4) [H3] runs in 0(n^ ) number of 
steps. Step (3) and (5) [HI] is done in 0(n) times. Step (6)[H4] runs in 
0(n^ ) number of steps Step (7) runs in 0(I*J) number of times and hence 
the whole heuristic mns in 0(n^ ) number of steps 
No. of iterations ; At each step of the heuristic we are assured of a 
strictly positive improvement in values of Wj’s : j gWJ (that can be raised 
in values ) and hence solution to problem DP can be improved by a 
positive quantity which is proportional to a quantity equal to 
min (Si, Sj) where, 

Ji = min(C22k-C2k*) 

k 

62 = niin (Extent_of_Decrease i) { as defined in 3.2.7 [H4] } 

1 


Let this number be denoted by S . And optimal solution to problem DP is 
a finite real number bound by sum of all C 1, j’s and C2j k’s . Hence 
maximum number of iterations required by heuristic H5 is given by [(sum 
of all Cl.j’s and C2 jl’s ) / ^] which is equal to c. 

Thus heuristic H5 runs in 0( c n^ ) number of times. 



CHAPTER 4 


COMPUTATIONAL EXPERIENCE AND ANALYSIS 

4.1 Introduction 

We have tested our heuristic on both small-size as well as large-size 
transshipment problems and in this chapter we give the computational results obtained by 
us from testing . 

As far as the computational complexity of the heuristic is concerned , we 
have already proved in the chapter three alongwith the heuristic solution procedure that 
our heuristic solves the transshipment problem in O (c n^ ) steps . Thus it has already been 
proved that our heuristic requires less computational effort as compared to the other 
available algorithms to solve the transshipment problem . Therefore in this chapter we put 
emphasis on the percentage closeness i e how close is the solution given by our heuristic 
to the optimal solution and the no. of iterations taken by our heuristic as compared by the 
linear programming optimization package CPLEX . Here we do not report the CPU time 
as a measure of our heuristic performance because CPU time depends greatly on finer 
details of the computational environment and the test problems such as : 

• Programming language , compiler , and computer ; 

• Implementation style and skill of the programmer ; 

• Generators used to generate the random test problems ; 

• Particular programming environment . 

Because of the multiple sources of variabilities , CPU times are often 
difficult to replicate , which is contrary to the spirit of scientific investigation . Another 


drawback of the use of CPU time is that it is an aggregate measure of empirical 
performance and does not provide much insight about an algorithm’s behaviour . For 
example , an algorithm generally performs some fundamental operations repeatedly, and a 
typical CPU time analysis does not help us to identify these ‘bottleneck’ operations 
Identifying the bottleneck operations of an algorithm can provide useful guidelines for 
where to provide future efforts to understand and subsequently improve an algorithm . 

All the transshipment problems we selected for testing are randomly 

generated . 

CPLEX Optimization Package 

For comparing the solution given by our heuristic with the optimal solution 
we use CPLEX Optimization Package available on system ‘BHASKAR’ CPLEX can 
take input in any of the format such as LP format , MPS format and SAV format . We 
used MPS input format and hence the input format of our heuristic is required to be 
converted into MPS format ( i.e column oriented format ) . During the solution process 
CPLEX reports its progress . The sblution process involves two stages . During Phase I, 
CPLEX searches for the feasible solution . In Phase n, CPLEX searches for the optimal 
feasible solution . The iteration log periodically displays the current iteration number and 
either the current scaled infeasibility during Phase I, or the objective function value during 
Phase II . This information can be useful for monitoring the rate of solution progress . The 
iteration log display can be modified using the ‘SET’ command to display differing 
amounts of data while the problem is being solved . Once the optimal solution has been 
found , the objective function value , solution time and iteration count are displayed by 


CPLEX . 


Computer Programs 


For testing our heuristic on different transportation problems, we prepare 
three different computer programs . First one is for generating random problems of 
specified sizes ( size ranging fi-om lOx 10x10 to 200 x 200 x 200 ) . The second one is 
for converting the input format for our heuristic into MPS format required by CPLEX. 
And the last one is for performing the step by step heuristic procedure . This last one is 
the main program of testing the solution procedure given by us in this thesis and it is given 
in appendix II 

Now in the following pages we give the results obtained by us in testing 
the various size problems as well as their percentage closeness to the optimal solution 
We start by giving some small sized problems followed by some medium size problems 
and finally we give the large size problems tested by us 


4.2 Small Size Transshipment Problems : 

We tested a total of 50 small sized problems The problem size in this 
category of problems is less than 10X10x10 ranging from 3x3x3 to 7x7x7. The 
distribution of the number of problems in this series of various sizes are as follows' 


Problem Size 

No. of Problems Taken 

Costs Range 

3x3x3 

15 

1 - 10* 

4x4x4 

15 

1 - 10* 

5x5x5 

10 

1 - 10* 

7x7x7 

10 

1 - 10* 


The result obtained by us in testing these small sized problems are given in 
Table 4.1 - Table 4.4 . 

The results coming from these testing show that the solution coming from 
these small size problems is optimal ( in most of the cases ) when compared with the 


results of CPLEX . 














TABLE 4.1 


PROBLEM SIZE : 3 x 3 x 3 


S. No. 

Solution Given by our 
Heuristic Solution 
Procedure 

Optimal Solution 
Using CPLEX 
Optimization Package 

Percentage 

Closeness 


Solution 

No. of 
Iterations 

Solution 

' 

No. of 
Iterations 


1 

11910.6 

6 

11911.356 

8 

99.99 

2 

27557.5 

7 

27557.5 

10 

100 

3 

16685.1 

6 

16685.1 

9 

100 

4 

22411.5 

6 

22411.5 

10 

100 

5 

20602.2 

7 

20602.2 

10 

100 

6 

22359.9 

6 

22359.9 

9 

100 

7 

25501.2 

7 

25501.2 

8 

100 

8 

23230.4 

8 

23230.4 

9 

100 

9 

23217.2 

7 

23217.2 

10 

100 

10 

23576.1 

7 

23576.1 

8 

100 

11 

15680.5 

7 

15680.5 

9 

100 

12 

16505.9 

8 

16505.9 

8 

100 

13 

12595.6 

5 

12595.6 

8 

100 

14 

27771.9 

7 

27771.9 

8 

100 

15 

21080.3 

6 

21080.3 

9 

100 





























































































TABLE 4.2 


PROBLEM SIZE ; 4 x 4 x 4 


S. No. Solution Given by our Optimal Solution Percentage 

Heuristic Solution Using CPLEX Closeness 

Procedure Optimization Package 

Solution I No. of Solution \ No. of 


Iterations Iterations 


1 

18105.5 

8 

18154.2 

12 

99.73 

2 

18519.8 

8 

18525.1 

12 

99.97 

3 

18526.0 

9 ' 


11 

100 

4 

16844.5 

8 

16844.5 

11 

100 

5 

13565.3 

8 

13565.3 

11 

100 

6 

15444.2 

8 

15444.2 

14 

100 

7 

18518.4 

10 

18518.4 

13 

100 

8 

82259.6 

13 

82259.6 

11 

100 

9 

10723.7 

10 

10723.7 

12 

100 

10 

10866.9 

8 

10866.9 

11 

100 

11 

15332.4 

10 

15332.4 

12 

100 

12 

83588.6 

10 

83588.6 

11 

100 

13 

11768.6 

9 

11768.6 

13 

100 

14 

12882.0 

8 

12882.0 

11 

100 

15 

17309.5 

9 

17309.5 

14 

100 


































































































TABLE 4.3 


PROBLEM SIZE : 5x5x5 


S. No. 

Solution Given by our 
Heuristic Solution 
Procedure 

Optimal Solution 
Using CPLEX 
Optimization Package 

Percentage 

Closeness 


Solution 

No. of 
Iterations 

Solution 

No. of 
Iterations 


1 

10096.6 

10 

10096.6 

16 

100 

2 

12176.2 

11 

12389.0 

17 

98.28 

3 

74788.2 

15 

7478.82 

16 

100 

4 

7670.39 

10 

7733.0 

15 

99.19 

5 

10290.6 

11 , 

10300.0 

14 

99.91 

6 

13214.0 

15 

13213.9 

15 

100 

7 

9401.56 

10 

9443.0 ! 

14 

99.56 

8 

11541.4 

11 

11541.43 

14 

100 

9 

8977.99 

11 

9118.0 

15 

98.46 

10 

12349.6 

16 

12355.0 

16 

99.96 



























































TABLE 4.4 


PROBLEM SIZE : 7x7x7 


S. No. 

Solution Given by our 

Optimal Solution 


Heuristic Solution 

Using CPLEX 


Procedure Optimization Package 



Solution 

No. of 
Iterations 

Solution 

No. of 
Iterations 

15656.2 

19 

15656.2 

27 

11766.0 

21 

11766.0 

24 

7671.85 

14 

7723.0 

23 

8434.69 

26 

8443.0 

25 

9247.58^ 

17 

9436.0 ^ 

25 

6207.0 

19 1 

6225.0 

22 

10054.4 

21 

10054.42 

20 

7351.84 

18 

7368.1 

21 

13356.0 

24 

13376.0 

23 

7079.91 

20 

7079.9 

20 


Percentage 

Closeness 


100 
100 
99.34 
99.90 
98.00 
99.71 
100 
99.78 ^ 

99.85 
100 
























































4.3 Medium Size Transshipment Problems 

In this category of transshipment problems , we randomly generate 20 
problems with size ranging from 10 x 10 x 10 to 75 x 75 x 75 . All these problems are 
generated using a random generator . The solution of these problems from our heuristic as 
well as that given by CPLEX are summarized in Table 4.5 . 


The problem size distribution in this category of problems is as follows • 


Problem Size 

No. Of Problems Tested 

Costs Range 

10 X 10 X 10 

3 

1 - 10* 

20 X 20 X 20 

3 

1 - 10* 

25 X 25 X 25 

3 

1 - 10* 

40 X 40 X 40 

3 

1 - 10* 

50x50x50 

3 

1 - 10* 

60 X 60 X 60 

2 

1 - 10* 

75 X 75 X 75 



3 

1 - 10* 


The results coming from these medium sized problems show that for this 
category of problems our solution is coming within 97-98 % to the optimal solution . 









TABLE 4.5 


PROBLEM OF SIZES RANGING FROM 10 x 10 x 10 to 75 x 75 x 75 


S. No. 


Solution Given by our 
Heuristic 


Optimal Solution Using 
CPLEX Optimization 
Package 



Solution 


6685.46 


6779.41 


5366.81 


No of 
Iterations 



1168.33 


1001.96 


1093.08 




Solution 


6927.67 


6961.3 


5406.013 


No. of 
Iterations 


32 


32 


29 


Percentage 

Closeness 



96.50 


97.38 


99.27 


4 

3931.32 

59 

4252.513 

70 

92.44 

5 


74 

4080.194 

73 

98.30 

6 

4082.95 

57 

4175.491 

69 

97.78 


7 

3341.71 

61 

3533.188 

90 

94.58 

8 

2687.67 

105 . 

2765.866 

84 

97.17 

9 

2385.34 

81 

2430.078 

83 

98.15 


10 

2101.39 

124 

2141.465 

140 

98.12 

11 

2085.19 

100 

2139.297 

131 

97.47 

12 

1840.84 

88 

1897.633 

149 

97.00 


13 

1470.34 

112 

1596.834 

168 

92.07 

14 

1663.51 

186 

1699.951 

192 

97.85 

15 

1628.35 

133 

1706.837 

175 

95.40 


16 

1423.45 

193 

1472.38 

208 

96.67 

17 

1322.62 

321 

1390.738 

206 

95.10 



97.10 


94.66 


95.94 














































































































4.4 Large Size Transshipment Problem 

In this category of transshipment problems also we generated 1 50 problem 
using random generator with size ranging from 100 x 100 x 100 to 200 x 200 x 200 . 
The solution obtained from our heuristic solution procedure as well as that given by 
CPLEX are summarized in Table 4.6 - Table 4.8 . Each table is for a particular size of 
transshipment problem and contains the summarized details of 50 problems . In each of 
these table we reported the problem size , Iterations taken by our heuristic and that of 
CPLEX and percentage closeness to the optimal solution for the solutions obtained using 
our heuristic solution procedure. The distribution of the various problems tested in this 
category is as follows ; 


Problem Size 

No. of Problems Taken 

Costs Range 

100 X 100 X 100 , 

50 

1-10* 

150 X 150 X 150 

50 

10 - 10000 

200 X 200 X 200 

50 

10 - 10000 


The results of the problems as summarized in the following tables show 
that our heuristic solution procedure for these large size problems gives a solution which 
is within 98 % (on average ) of the optimal . 















TABLE 4.6 


PROBLEM SIZE : 100x100x100 

S. No. No. of Iteration No. of Iterations Percentage 
Taken by our Taken by CPLEX Closeness to the 
Heuristic to give optimization Package Optimal Solution 
the Solution to give the Optimal by our Heuristic 

Solution 



96.57959 


96.37769 


95.20112 


97.25020 


96.81606 


96.25575 


96.23152 


98.32666 


96.01851 


96.73304 


101.8670 


95.65126 


97.36515 


95.53120 


96.29604 


96.75471 


99.77138 


95.52540 


94.91795 


96.40491 


98.17961 


96.98930 


98.78578 


94.69378 


94.69325 
























































































TABLE 4.6 [Contd ...] 


S. No. No. of Iteration No. of Iterations Percentage 

Taken by our Taken by CPLEX Closeness to the 
Heuristic to give optimization Package Optimal Solution 
the Solution to give the Optimal by our Heuristic 

Solution 


26 

262 

27 

349 

28 

336 

29 

370 

30 

222 

31 

213 

32 

329 

33 

268 

34 

283 

35 

318 

36 

329 

37 

310 

38 

303 

39 

315 

40 

282 

: 41 

415 

1 42 ' 

350 

43 

336 

44 

370 

45 

283 

46 

291 

47 

281 

48 

308 

49[aal] 

435 

50 

285 



Average No. of Iterations Taken by Our Heuristic 
Average No. of Iterations Taken by CPLEX 
Average Percentage Closeness to the Optimal Solution 


96.98930 


98.78578 


94.69325 


97.58082 


95.25702 


96.39582 


95.52060 


96.98929 


94.69325 


96.98982 


97.02474 


94.69325 


95.82753 


95.82742 


97.58093 


98.17963 


94.69325 


97.58082 


95.25645 


97.58083 


96.40491 


97.58083 


96.39582 


96.98932 


95.86452 


316 
377 
96.6 % 
































































































TABLE 4.7 


PROBLEM SIZE : 150 x 150 x 150 


S. No. 

No. of Iteration 
Taken by our 
Heuristic to give the 
Solution 

No. of Iterations Taken 
by CPLEX 
optimization Package 
to give the Optimal 
Solution 



Percentage 
Closeness to 
the Optimal 
Solution by 
our Heuristic 

1 

339 

601 

98.14 

2 

328 

419 

93.05 

3 

396 

635 

98.79 

4 

356 

598 

98.61 

5 

352 

588 

98.71 

6 

388 

737 

99.29 

7 

325 

595 

98.80 

8 

337 

628 

99.14 

9 

364 

895 

100.27 

10 

367 

591 

98.44 

11 

354 

806 

98.89 

12 

393 

472 

98.87 

13 

373 

472 

98.87 

14 i 

395 

659 

97.07 

15 1 

380 

697 

97.24 

16 ! 

370 

509 

98.84 

17 

353 

547 

97.08 

18 

359 

472 

96.72 

19 

377 

584 

97.73 

20 

380 

549 

95.90 

21 

371 

734 

96.68 

22 

388 

513 

97.89 

23 

354 

731 

97 99 

24 

393 

622 

97.58 

25 

384 

685 

97.43 


























































































TABLE 4.7 [Contd....] 


S. No. 

No. of Iteration Taken 
by our Heuristic to 
give the Solution 

No. of Iterations Taken by 
CPLEX optimization 
Package to give the 
Optimal Solution 

Percentage 
Closeness to the 
Optimal Solution 
by our Heuristic 

26 

351 

697 

96.72 

27 

377 

734 

97.78 

28 

365 

584 

96.37 

29 

366 

622 

97.43 

30 

381 

584 

96.37 

31 

395 

510 

97.07 

32 

358 

735 

97.78 

33 

370 

697 

96.72 

34 

388 

659 

98.51 

35 

377 

734 

98.47 

36 

383 

584 

96.66 

37 

392 

472 

97.28 

38 

373 

713 

95.73 

39 

364 

553 

97.59 

40 

381 

665 

97.48 

41 

370 

622 

97.07 

42 

383 

659 

99.42 

43 

370 

473 

96.05 

44 

393 

549 

92.17 

45 

359 

697 

99.18 

46 

380 

673 

96.72 

47 

377 

624 

97.43 

48 

381 1 

697 

96.72 

49 

359 

472 

98.87 

50 

383 

628 

98.81 


Average No. of Iterations Taken by our Heuristic 
Average No. of Iterations Taken by CPLEX 
Average Percentage Closeness to the Optimal 


371 
620 
97.5 % 


















































































TABLE 4.8 


PROBLEM SIZE : 200x200x200 


S. No. No. of Iteration No. of Iterations Taken Percentage 
Taken by our by CPLEX optimization Closeness to the 
Heuristic to give Package to give the Optimal Solution 
the Solution Optimal Solution by our Heuristic 



1458 


1702 


1601 


1415 


1741 


1645 


1299 


1388 


1846 


1371 


1705 


1366 


1820 


1625 


1690 


1495 


1755 


1560 


1820 


1690 


1625 


1560 


1495 


93 


99.79 


100.60 


99.56 


99.40 


100.00 


99.88 


100.26 


99.59 


100.12 


99.96 


99.43 


99.09 


98.63 


99.46 


100.15 


100.21 


99.65 


100.13 


100.74 


101.50 


100.36 


99.24 


102.03 


99.33 


99.67 



















































































TABLE 4.8 [ Contd ....] 


S. No. 

No. of Iteration 
Taken by our 
Heuristic to give 
the Solution 

No. of Iterations Taken 
by CPLEX optimization 
Package to give the 
Optima! Solution 

Percentage 
Closeness to the 
Optimal Solution 
by our Heuristic 

26 

673 

1820 

100.07 

27 

873 

1625 

99.47 

28 

757 

1823 

98.95 

29 

674 

1365 

99.18 

30 

555 

1560 

99.24 

31 

878 

1755 

99.61 

32 

491 

1629 

99.46 

33 

672 


101.04 

34 

872 1 

1365 

99.15 

35 

708 

1820 

98.95 

36 1 

508 

1560 

99.24 

37 

735 

1690 

101.50 

38 

645 

1893 

100.07 

39 

809 

1690 

98.81 

40 

709 

1560 

99.24 

41 

800 

1625 

100.36 

42 

799 

1690 

100.60 

43 

491 

1365 

99.16 

44 

717 

1625 

100.36 

45 

682 

1430 

99.98 

46 

532 

1820 

100.74 

47 

626 

1627 

100.36 

48 

757 

1690 

101.51 

49 

631 

1495 

99.31 

50 

777 

1389 

99.83 


Average No. of Iterations Taken by our Heuristic 
Average No. of Iterations Taken by CPLEX 
Average Percentage Closeness to the Optimal Result 


680 
1604 
99.3 % 






















4.5 A Note on the No. Of Iterations Taken by Our Heuristic i 

As seen from the Tables given in the previous pages as well as the charts 
attached at suitable places we can very well see that the number of iterations taken by our 
solution procedure is considerably low than that was obtained from the CPLEX even for the 
large size problems . 

Though for the small and medium size problem the number of iterations coming 
from our solution procedure is slightly lower than that of the CPLEX , the clear picture in the 
difference in number of iterations taken can be seen in large size problems 

• For problem size 100 x 100 x 100 , it takes iterations 0.75-0.80 (on average) times the 
iterations taken by the CPLEX ; 

• For problem size 150 x 150 x 150 , it takes iterations 0.50-0.60 (on average) times the 
iterations taken by the CPLEX ; 

• For problem size 200 x 200 x 200 , it takes iterations 0.40-0.50 (on average) times the 
iterations taken by the CPLEX ; 

Thus we can see that our solution heuristic procedure come to be efficient in its 
solution as well as in the number of iterations taken by it when compared with the CPLEX 


optimizing package . 



4.6 Running our Solution Procedure in 0( c n^) number of steps : 

Though our solution procedure runs in 0(c n^) number of steps, now we 
show the results coming from our heuristic procedure if we run it in 0(c n^) number of 
steps As already explained in chapter three all but Heuristic H3 [3.5.6] runs in either 0(n) 
times or in O(n^) times . Only Heuristic H3 [3.5.6] runs in 0( n^) times . So if we 
eliminate the improvement in our solution coming from heuristic H3 , i.e if we don’t 
include improvements by H3 , we can successfully run our heuristic procedure in 0( c n^ ) 
number of times. 

We solved 150 problems of sizes varying from 100 x 100 x 100 to 200 x 
200 x 200 to show the result coming from running our solution procedure both in 0(c n^) 
and 0( c n^) times In the Tables 4.9 - 4.1 1 , percentage closeness to the optimal solution 
both by running the solution procedure in 0(c n^) and 0( c n^) times are shown 

The result shows that even in runniag the heuristic in 0(c n^) times, we get 
good results which are also very close to the optimal solution. 



TABLE 4.9 

Problem Size : 100 x 100 x 100 


Solution After H2 only 
Running in O(cn^) 


Iterations 

Percentage 

Closeness 

248 

97.22 

335 

96.43 

311 

96.64 

263 

95.98 

298 

96.92 

284 

97.05 

224 

96.99 

248 

96.03 

298 

97.09 

347 

96.67 

272 

96.20 

283 

96.43 

218 

96.61 

227 

95.98 

221 

96.77 

251 

96.79 

326 

96.91 

223 

96.66 

224 

96.45 

371 

96.36 

400 

96.44 

187 

96.68 

197 

96.69 

187 

96.89 

268 

96.81 


Solution After H2+H3 
Running in O(cn^) 

Iterations 
taken by 
CPLEX 

iterations 

Percentage 

Closeness 


307 

97.31 

535 

353 

96.47 

534 

337 

96.68 

380 

303 

96.08 

400 

329 

96.96 

377 

289 

97.07 

390 

264 

97.05 

431 

268 

96.06 

419 

315 

97.11 

586 

373 

96.71 

401 

307 

96.26 

423 

286 

96.43 

365 

269 

96.66 

367 

291 1 

96.10 

422 

277 

96.88 

468 

302 

96.90 

613 

341 

96.94 

531 

274 

96.72 

380 

323 

96.66 

435 

425 

96.44 

406 

463 

96.52 

521 

227 

96.76 

601 

216 

96.72 

429 

204 

96.91 

484 

297 

96.86 

399 


: ’"/oa 
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TABLE 4.9 [Contd....] 


Solution After H2 only 
Running in O(cn^) 

Solution After H2+H3 
Running in O(cn^) 

Iterations 

Percentage 

Closeness 

Iterations 

Percentage 

Closeness 

275 

96.44 

342 

96.52 

305 

96.15 

333 

96.20 

318 

96.63 

331 

96.67 

287 

96.11 

298 

96.14 

242 

96.31 

277 

96.35 

267 

97.06 

288 

97.09 

233 

96.43 

253 

96.46 

276 

96.68 

300 

96.73 

262 

96.89 

307 

96.96 

251 

96.56 

293 

96.67 

315 

95.29 

331 

95.30 

263 

97.05 

285 

97.09 

238 

96.80 

279 

96.86 

241 

96.80 

265 

96.83 

246 

96.48 

283 

96.55 

288 

96.78 

310 

96.81 

310 

97.26 

354 

97.32 

236 

97.07 

317 

97.17 

261 

96.85 

319 

97.00 

357 

96.55 

383 

96.59 

291 

96.41 

342 

96.51 

186 

96.94 

263 

97.06 

237 

96.65 

278 

96.72 

258 

96.60 

304 

96.69 

274 

96.53 

293 

96.57 


Iterations 
taken by 
CPLEX 



384 


407 


454 


376 


413 


589 


378 


374 


392 


400 


373 


456 


574 


418 


365 


400 


497 


617 


398 


503 


397 


380 


378 


559 


363 
































































































































TABLE 4.10 

Problem Size : 150 x 150 x 150 


Solution After H2 only I Solution After H2+H3 
Running in ©(cn^) Running in O(cn^) 


Iterations 
taken by 
CPLEX 


iterations 



5 

339 

6 

371 

7 

348 

8 

405 

9 

323 

10 

372 

11 

352 

12 

400 

13 

386 

14 

366 

15 

343 

16 

362 

17 

313 

18 

328 

19 

380 

20 

384 

21 

351 

22 

375 

23 

343 

24 

384 

25 

328 


Percentage 

Closeness 


98.43 


.98.39 


99.18 


98.10 


98.35 


97.92 


99.03 


99.14 


99.31 


98.20 


99.45 


99.67 


99.46 


99.55 


98.63 


98.56 


99.55 


95.39 


98.52 



iterations 

Percentage 

Closeness 

360 

98.45 

371 

98.51 

341 

99.23 

378 

98.35 

345 

98.45 

383 

98.45 

353 

99.08 

410 

99.21 

342 

99.51 

383 

99.21 

369 

99.73 

406 

99.73 

388 

99.51 

376 

99.67 

357 

98.82 

392 

98.96 

320 

99.61 

351 

96.04 

380 

98.52 

392 

99.50 

351 

99.06 

384 

99.16 

343 

99.21 

392 

99.40 

340 

98.75 



623 


604 


626 


565 


582 


613 


565 


640 


716 


578 


618 


561 


641 


575 


636 


675 








































































































































TABLE 4.10 [Contd...] 


S.N. 


Solution After H2 only 
Running in O(cn^) 



iterations 


309 


355 


381 


346 


347 


372 


370 


380 


383 


362 


345 


379 


376 


375 


362 


364 


326 


356 


346 


370 


Percentage 

Closeness 


97.88 


98.89 


. 99.58 


99.46 


99.09 


98.86 


99.21 


98.32 


98.80 


99.71 


98.98 


98.95 


99.77 


99.09 


99.03 


98.95 


98.76 


98.94 


99.46 


99.09 


99.18 


96.79 


98.50 


95.63 


99.57 


Solution Atter H2+H3 It 
Running in 0(cn^) t 

iterations 

Percentage 

Closeness 

325 

98.03 

379 

99.34 

381 

99.58 

362 

99.69 

364 

99.33 

377 

98.96 

378 

99.30 

401 

98.59 

383 

98.80 

362 

99.71 

369 

99.38 

399 

99.12 

388 

99.92 

382 

99.18 

366 

99.17 

370 

99.04 

354 

99.28 

356 

98.94 

362 

99.69 

375 

99.20 

384 

99.41 

326 

97.59 

383 

98.72 

386 

97.36 

393 

99.63 


CPLEX 


595 


535 


573 


587 


595 


635 


658 


604 


707 


624 


582 


644 


601 


577 


610 


613 


619 


606 


587 


567 


593 


595 


589 


696 


597 



























































































































TABLE 4.11 

Problem Size : 200 x 200 x 200 


S.N. 


Solution After H2 only 
Running in O(cn^) 


Solution After H2+H3 
Running in O(cn^) 


Iterations 
taken by 
CPLEX 



Iterations 


289 


711 


283 


521 


623 


457 


711 


758 


480 


410 


815 


270 


494 


446 


503 


917 


694 


275 


930 


452 


569 


535 


364 


312 


375 


Percentage 

Closeness 


99.98 


99.19 


99.97 


99.96 


99.86 


99.76 


99.88 


99.88 


99.19 


99.88 


99.68 


99.96 


99.16 


99.96 


99.53 


99.67 


99.65 


Iterations 


289 


711 


283 


521 


623 


457 


711 


758 


480 


410 


815 


270 


494 


446 


503 


917 


694 


275 


930 


452 


569 


535 


364 


312 


375 


Percentage 

Closeness 


99.98 


99.19 


99.97 


99.96 


99.86 


99.76 


99.88 


99.88 


99.19 


99.88 


99.68 


99.96 


99.16 


99.96 


99.53 


99.67 


99.65 


99.98 


99.29 


99.65 


99.85 


99.66 


99.98 


99.78 


99.52 



1281 


1599 


1620 


1850 


1762 


1569 


1497 


1260 


1315 


1554 


1596 


1463 


1159 


1315 


1346 


1441 


1340 


1286 


1780 


1410 


1603 


1328 


1320 


1618 


1330 






































































































TABLE 4.11 [Contd...] 


S.N. 


Solution After H2 only 
Running in O(cn^) 


Iterations 


Percentage 

Closeness 


26 


787 


99.98 


27 


478 


99.56 


28 


476 


99.66 


29 


502 


99.67 


30 


417 


99.98 


31 


345 


99.99 


32 


299 


99.48 


33 


361 


99.54 


34 


539 


99.64 


35 


628 


99.96 


36 


890 


99.55 


37 


201 


99.65 


38 


284 


99.81 


39 


439 


99.39 


40 


201 


99.65 


41 


495 


99.66 


42 


616 


99.67 


43 


845 


99.66 


44 


628 


99.64 


45 


628 


99.66 


46 


283 


99.61 


47 


585 


99.76 


48 


591 


99.55 


49 


50 


486 


583 



Solution After H2+H3 
Running in O(cn^) 

iterations 
taken by 
CPLEX 


Iterations 

Percentage 

Closeness 



787 

99.98 

1266 


478 

99.56 

1432 


476 

99.66 

1413 


502 

99.67 

1307 


417 

99.98 

1803 


345 

99.99 

1830 


299 

99.48 

1817 


361 

99.54 

1619 


539 

99.64 

1462 


628 

99.96 

1424 


890 

99.55 

1266 


201 

99.65 

1281 


284 

99.81 

1442 


439 

99.39 

1263 


201 

99.65 

1466 


495 

99.66 

1370 


616 

99.67 

1245 


845 

99.66 

1388 


628 

99.64 

1354 


628 

99.66 

1366 


283 

99.61 

1361 


585 

99.76 

1344 I 


591 

99.55 

1241 


486 

99.57 

1571 


583 

99.51 

1381 


99.57 

99.51 




























































































































CHAPTER 5 


CONCLUSIONS 

AND 

DIRECTIONS FOR FUTURE RESEARCH 

As far as the computational complexity of the solution procedure given by us is 
concerned , we can say that our heuristic runs in 0(c n^) number of step . Using this solution 
procedure we solved transshipment problems of small as well as large size . It is found that 
even for the large sized problems the solution comes closer to 97 % to the optimal solution 
obtained by CPLEX . We also showed that even while running the heuristic in 0( c n^) times , 
the solution procedure gives good results . Thus the solution coming from our solution 
procedure can give a good start to conventional dual based algorithms to reach the optimal 

solution. 

We hope that future work would extend our work to multistage minimum cost 
flow problems with no restrictions on capacity of flow on arcs. May be later work can 
concentrate on capacitated network flows problem as well 
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AmLrNuiA 1 


Some Solved Numerical Problems Using the Heuristic 

In order to fully illustrate the working of the solution procedure given by us in this 
paper , we are here giving three sample problems . 

Problem 1 : 

Problem Size : 3 X 3 X 3 
Given ; 

Supplies b, =. (0.2 , 0.5 , 0.3) 

Demands dk = (0 25 , 0 4 , 0.35) 

■3 5 .8 

Matrix Cl ,j= 9 6 1 

1 3 8 

'2 1 4 

Matrix C2jk= 3 2 5 

6 7 6 


C2ji = (2,3,6) 
C2 j2=(i;2,7) 
C2 j3 = (4,5,6) 
C2k* = (2, 1 ,4) 

12 ={ 1 } 

13 = {1} 

V2o = 4 


Solution Procedure : 

Iteration # 0 ; 

Cl ij = a,5,8) 
Cl 2j = (9,6,i) 
Cl 3 j = (2,3,8) 
Cir = (3,l,2) 

mi ={1} 
m2={3} 
m3={l} 

Thus we have , 

Vlo = 3 


61 



Zi = 0, Z 2 = 2, Z 3 = 1 

Vi=2, V2 = 3, V3 = 0 

Zdual = (Vlo - ^ Z, b, ) 
1=1 

+ (V2o- i Vk.dk) 

k=l 

= 4 00 


Iteration # 1 


WJ = 1, = {1} 

using H2 

gain = di + d 2 + ds = 1.0 

loss = bi + bs = 0.5 

^®t_Benefit — gain - loss — 0.5 [>0 Therefore Increment Possible ] 

Extent_of_Increase2 = min ( 1 , 1 , 1) = 1 

Extent_of^decreasel = min ( 

00 , 8 , 00 ) = 8 

Extnet_of^Increase - min (Extent_Oi^Increase2, Extent_oy)ecrease2) = 1 

Cl ij = a,5,8) 

C2ji=(3,3,6) 

Cl2j = ( 8 , 6 ,i) 

C2j2=(2,2,7) 

Cl3j = (i,3,8) 

C2j3=(5,5,6) 

Cli* = (2, 1 , 1) 

C2k* = (3,2,5) 

mi ={1} 

li = {l,2} 

m2={3} 

12={1,2} 

m3={l} 

13 = {1,2}. 

Thus we have , 


Vlo = 2 

V2o = 5 

Z] = 0, Z 2 = 1, Z 3 = 1 

Vi=2, V2 = 3, V3 = 0 

Zdual =(Vlo-Xz, b.) 

1=1 

+ (V2o - i: Vk.dk) 

k=l 

= 4.50 


Iteration # 2 


WJ = Ii = { 1,2} 

using H2 _ 

gain = di+d2+d3 = l.C 

1 loss = bi+bs = 0.5 

NetJBenefit = gtdn - 1 

loss = 05 [>0 Therefore Increment Possible ] 


an 



Extent_of_Increase2 = min ( 3 , 5 , 1) = 1 
Extent_of_decreasel = min ( oo,5, oo) = 5 

Extnet_ofJncrease = min (Extent_of_Increase2, Extent_of_Decrease2) = 1 


Cl ij = a,4,8) 
Cl 2j = (7,5.i) 
Cl3j = (g,2,8) 

Cli‘ = (l,l,0) 

mi ={1} 
m2={3} 
m3={l} 

Thus we have , 

Vlo= 1 

Zi = 0, Z2 = 0, Z3 = 1 


C2ji=(4,4,6) 
C2j2=(3,3,7) 
C2j 3=(6,6,6) 
C2; = (4,3,6) 
li = {l,2} 

12 ={ 1 , 2 } 
l3 = {1,2,3} 

V2o = 6 

Vi=2, V2 = 3, V3 = 0 


^ K 

Z dual ~ (Vlo ~ ^ z, . bi ) + (V2o - 5^ Vfc . dt ) 

1=1 k=l 

= 5.00 

Iteration # 3 : 

WJ = li = {l,2} using H2 
gain = di+d 2 = 0 65 loss = bi+b 3 = 0.5 

Net_Benefit = gain - loss = 0. 15 [>0 Therefore Increment Possible ] 
Extent_of_Increase2 = min ( 2,3, oo) = 2 
Extent_of_decreasel = min ( oo,4, oo) = 4 

Extnet_o^Increase = min (Extent_of_Increase2, Extent_of_Decrease2) = 2 


Cl ij = (-l,2,8) 
Cl 2j = (5,3,i) 
Cl3j = (=2,0,8) 
Cli* = (-l,l,-2) 

mi ={1} 


C2 


6, 

6) 

C2 

32=(5^ 

5, 

7) 

C2 

j3=(8. 

8, 

6) 

C2 

»:* = (6, 

5,' 

5) 

li = 

= {1,2,3} 




✓'--I 


m2={3} 

m3={l} 


12 ={ 1 , 2 } 

I 3 = {3} 


Thus we have , 

Vlo =1 V2o = 6 

Zi = 2, Z 2 = 0, Z 3 = 3 Vi = 0 , V 2 = 1 , V 3 = 0 

Z dual ~ (Vlo - Z, . b, ) + (V2o - 5] Vk.dk) 

1=1 k=l 

= 5.30 

Iteration # 4 

WJ = li = { 1,2 } using H 2 

gain = di+d 2 +d 3 = 1 .00 loss = bi+b 2 +b 3 = 1.00 

Net_Benefit = gain - loss = 0 [ Therefore No Increment Possible ] 

Other combinations of the current Ik’s also using H 2 shows Net_Benefit <= 0 , 
hence now we go for finding WJ using H3. 

Iteration # 5 ; 

WJ = {1,2,3} ' using H3 

gain = loss = 1 0 Hence no improvement possible 
Other combinations of Ik’s also show Net_benefit<=0 , hence no improvement 
possible at all. 

Heuristic Terminates Here 

Objective function Value Zduai = 5.30 
No. of iterations taken = 5 

Optimal Results From CPLEX = 5.30 


Problem 2 : 


A 


Problem Size : 4 X 4 X 4 
Given 


Matrix Cl ij= 


Matrix C2 j k • 


1 . 


1,0.3, 0.5, 0.1) 

2,0.3,04,0 1) 

. ...2.... 

. i 6 

3.. . 

...4... .2 

5.... 

...8 6 

6. .. 

...2 1 

3.... 

...7 8" 

9 .. 

...5... .3 

. . ..3. .. 

..2 . 4 

3 

..2 3 


Solution Procedure : 
Iteration # 0 ; 


Thus we have, 


Cl ,j = (3,2,5,6) 

C2j,= (2,6,l,4) 

Cl2j = (4,3,4,2) 

C2j2=a,9,3_,3) 

Cl 3j = (2,5,8,6) 

C2j3=(7,5,J^ 

ci4j = (7,6.,2, ri 

C2j4=(8,14,3) 

C1>(2,2,2,1) 

C2k* = (l,3,2,3) 

mi ={2} 

li = {3} 

m2 ={4} 

J2={1,3,4} 

m3={l} 

l3 = {3,4} 

m4={4) 

l4={2,4} 

Vlo = 2 

V2o = 3 

Z2 = 0, Z3 = 0, Z4 = 1 

Vi = 2, V2 =0,V3 =1,V4 = 0 


= (Vlo-i; Zi.b,) + (V2o- s Vk.4) 

k=l 

1 = 1 



Iteration # 1 


WJ - 1 1 - { 3 ) using H2 
gain = d, =0.2 loss = 0 

Net_Benefit = gain - loss = 02 [>0 Therefore Increment Possible ] 
Extent_of_Increase2 = min ( 1, oo, oo, oo) = 1 
Extent_of_decreasel = min ( 3 , 2, 6 , 1) =1 

Extnet_of_Increase = min (Extent_o^Increase2, Extent_o^Decrease2) = 1 


Cl ij = (3, 2,4,6) 

C2ji=(2,6,2,4) 

Cl2j = (4,3,3,2) 

C2j2=a,9,4,3) 

Cl3j = (2,5,7,6) 

C2,3=(7,5,3,^ 

Cl4j = (7,6,l_l) 

C2 j4=(8-,3,5,3) 

Cl.‘ = (2,2,2,l) 

C2k* = (2,3,2,3) 

mi ={2} 

li = {1.3) 

m2 ={4} 

12={1, 4} 

ms ={1} 

l3 = {4} 

m4={3 , 4} 

l4={2,4} 

11 

O 

> 

V2o = 3 

Z2 = 0, Zs = 0, Z4 = 1 

Vi = 1, V2 =0,V3 =1,V4 = 0 


r=l 


Zdua. = (Vlo-1; Zi.b.) + (V2o- E Vk.dt) 

k=l 


= 4 30 

Iteration # 2 : 

WJ = li = {l,3} usingrn 

gain = di = 0.2 loss = bs + bd = 0.6 

NetJBenefit = gain - loss = -0.4 [<0 Therefore No Increment Possible ] 
Other combinations of the current Ik’s also show Net_Benefit <=0 using H2 
and hence we now go to find WJ using H3. 

Iteration # 3 : 
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WJ = {1,3, 4} using H3. 

gain = di+d2+d3 = 0.9 loss = b 2 +b 3 +b 4 = 0.9 

Net_Benefit =0 Therfore no Improvement possible. 

Other combinations of the current It’s also show no improvement using H3. 
Hence no improvement at all 
Heuristic terminates here 

Objective function value Zduai = 4.30 
No of Iterarionstaken =3 
Optimal Result from CPLEX =4.30 


Problem 3 : 

Problem Size ; 3 X 4 X 3 
Given . 


Supplies b, = (0.3 , 0.5 , 0.2 ) 
Demands d t = (0 1 , 0.7 , 0.2 ) 


Matrix Cl 


3 .4 2 

5 6 4.. 

3 2 5. 


.3 

.1 

.3 


Matrix C2 j t = 


3 2 5 

6 1 3 

2 1 4 

4 3 1 


Solution Procedure : 
Iteration # 0 ; 


Cl 1, = (3,4, 2, 3) 

C2 ji = (3,6, 2,4) 

Cl2j = (5,6,4, i) 

C2 j2=(2, i,i,3) 

Cl3j = (3, 2,5,3) 

C2 j3 = (5,3, 4,1) 

Cir = (2,l,2) 

C2t* = (2,l,l) 

mi ={3) 

li = {3} 

m2={4} 

l2={2,3} 



ni3=i2) 


l3 = {4} 


Thus we have, 


Vlo = 2 

Zi = 0, Z 2 = 1 , Z 3 = 0 


V2o = 2 

Vi = 0, V2 =1,V3 =1 


z 


dual 


= (V1.-E 


2i-b,) + (V2o- I v,.dk) 

k=I 


= 2.60 

I teration #I 

WJ = 1, = {3} usingm 
gain = di = 0.1 loss = 63 = 0.3 

Net_Benefit = gain - loss = - 0.2 [<0 Therefore No Increment Possible ] 


Iteration # 2 : 


WJ = I 2 = { 2 , 3 } using H 2 . 
gain = di+d 2 = 0.8 loss = 61 + 63 = 0.5 

Net_Benefit = gain - loss = 0.3 [>0 Hence Improvement possible] 

Extent_of_Increase 2 = min (00, 1 , 00) = 1 
Extent_of_decreasel = min (00, 3, 00) = 3 

Extnet_of_Increase = min (Extent_ofJ[ncrease 2 , Extent_of_Decrease 2 ) 
Clij = (3,3,i,3). C2ji = (3,7, 3,4) 

Cl2j = (5,5,3, i) C2j2=(2, 2,2,3) 

Cl3j = (3, 1,4,3) C2j3 = (5,4, 5,1) 

Ci; = (l,l,l) C 2 ; = (3,2 ,1) 


mi ={3} 
m2={4} 
m 3 ={ 2 } 


12={1,2,3} 

l3={4} 


Thus we have, 

Vlo= 1 

zi = 0, Z2 = 0, Z3 = 0 


V2o = 3 

Vl = 0, V2 =1,V3 =2 
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Zduai -(Vlo-^ Zi •bi) + (V2o- 2 ; Vk.dk) 

i=I k=l 


= 2.90 

Iteration #3 , 

WJ = 1 1 = { 1,3} using H2 

gain = di =0.1 loss =bi = 0.3 

Net_Benefit = gain - loss = - 0.2 [<0 Therefore No Increment Possible ] 


Iteration # 4 


WJ = l2= {1,2,3} 

gcun = di+d 2 = 0.8 loss = bi+b 3 = 0.5 

Net_Benefit = gain - loss = 0 3 [ >0 Therefore Improvement possible] 

Extent_of_Increase2 = min ( 00 , 1 , 00 ) = 1 

Extent_of_Decreasel = min ( 00 , 2 , 00 ) = 2 

Extent_of_Increase = min (Extent_of_Increase2 , Extent_of_Decrease2) 


Cl ij = (2,2,0,3) 

Cl2j = (4,4,2, i) 

Cl3j = (2, 0,3,3) 

C1>(0,1,0) 

mi ={3} 

m2={4} 

m3 ={2} 

Thus we have, 

Vlo= 1 


C2ji=(4,8, 4,4) 
C2 j2=(3, 3,3,3) 
C2j3=(6,5, 6,1) 
C2k* = (4,3 ,1) 
li = {l,3,4} 
I2={1,2,3,4} 
l3 = {4} 

V2o = 4 


Zi = 1,Z2 = 0, Z3= 1 Vi = 0, Vi=l,Vi=3 


^ dual = (Vlo - X 2 , . b, ) + (V2o - z Vk.dk) 


= 3.20 

Iteration # 5 : 
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WJ = 1 1 = { 1,3,4} using H2 

gain = di + ds = 0.3 loss = bi + b 2 = 0.8 

Net_Benefit = gmn - loss = - 0.5 [<0 Therefore No Increment Possible ] 
Other combinations of the current Vs also using H2 shows NetJBenefit <=0 , 
hence now we go for finding WJ usihg H3. 

The combinations of the current Vs also show no improvement using H3. 
Hence no Improvement at all. 

Heuristic terminates here 

Objective function value Zduai = 3.20 
No. of Iterarions taken = 5 


Optimal Results from CPLEX = 3.20 



APPENDIX n 


Pascal Program made to test the solution procedure given in this paper for 
various randomly generated Problems . 

Program Heuristic(Input, Output), 

Label 100,200, (* Lables Declaration *) 

Const (* Constants Declaration *) 

XX = 200, 
yy = 200, 
zz = 200; 
inf= 100000000, 


Type 

Digits = set of 1 ..200, 


Var (* Variables Declaration *) 

over • boolean; 

dec,Cl . array[l. xx,l..yy] of integer; 

C2 : array[l..yy,l..zz] of integer, 
b : array[l. xx] of real, 
d : array[l..zz] of real; 
m : array[l xx] of digits; 

1 array[l zz] of digits, 

kl , w, : digits; 

u , i , j , k , k2 integer; 

cll , clstar : array[l..xx] of integer; 

c22 , c2star : array[l..zz] of integer; 

eoi : array[l..zz] of integer; ' 

eod ; array[l..xx] of integer; 

eoil , eodl , feoi • integer; 
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X , y , z : integer; 

opt 5 gain , loss , net_benefit real, 
iterations ; integer; 

5oG«Go()e(5o()e()c()(()e()eGe()oO(e<)oOcGo(5ooooo()ooo(3cGOOtOoOeOoOoOeGoO^^ 

Procedure To Initialize all the variables used in the Program 

Procedure Initialize, ' 

var 

ij,k:integer; 

Begin 

for i:=l to xx do 
for j:=l to yy do 
cl[ij];=0; 
forj =1 to yy do 
for k:=l to zz do 
c2[j,k];=0; 
for i:=l to xx do 
b[i]:=0.0; 
for k:=l to zz do 
d[k];=0.0; 
for i:=l to xx do 
m[i]:=[]; 
for k;=l to zz do 

I[k]=[], 

gain:=0.0; 

loss;=0.0; 

net_benefit:=0.0; 

iterations:=0; 

End; 

)0C)e()!O0()^)P?)0C)0()8G0O0OPO0OeO0O9OeO0O«e()cC)0GflOeO0O8O0()JG0Ot^^ 
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Procedure to Read Input for the Problem 


Procedure Read_input, 
var 

i , j , k ; integer. 

Begin 

read(x), (* Supply Nodes *) 

read(y); (* Intermediate Nodes *) 

read(z), (* Demand Nodes *) 

for i:=l to x do 

read(b[i]), (* Supplies *) 

for k:=l to z do 

read(d[k]), (* Demands *) 

for i:=l to x do 
begin 

for j:=l to y do 

read(cl[i,j]), (* Unit Transportation Cost *) 

readin, 

end; 

forj:=ltoydo 

begin 

for k:=l to z do 

read(c2[j,k]), (* Unit Transportation Cost *) 

readin, 

end; 

End; 


Procedure to Find the Values of Cl ,* and C2t’ 


Procedure Find_Stars, 
var 

ij,k: integer, 


Begin 

for i:=l to X do 
begin 

clstar[i].=cl[i,l]; » 
for j:=l to y do 
begin 

if (clstar[i]>cl[ij]) then 
clstar[i]:=cl[ij]; 

end, 

end; 

for k =1 to z do 
begin 

c2star[k] =c2[l,k]; 
for j;=l to y do 
begin 

if (c2star[k]>c2|j,k]) then 
c2star[k]:=c2|j,k], 

end, 

end; 

End; 

)oOOOe()(()o()o()OGe()eO!()O()8()«C)O()cOO()JGoGeOoOO0OO«OoOoG(OJ()OG«GcGoOOG«Of(5OC)!GOG«OoO(O«GC5CGflOflGOGOG<GO( 



Procedure to Find the Sets mi and I ^ used in the program 


Procedure find_sets_m_andj, 
var 

i , j , k . integer. 

Begin 

for i:=l to x do 
m[i]:=[]; 

for k;=l to z do 
l[k]:=[]; 

for i:=l to x do 

for j;=l to y do 

if (cl[ij]=clstar[i]) then 
m[i]:=m[i]+[j]; • 

for k:=l to z do 

for j:=l to y do 

if (c2[j,k]=c2star[k]) then 

l[k]:=l[k]+D]; 


End; 
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Procedure to Find Net Benefit 


Procedure Find_net_benefit; 
var 

ij , k ; integer, 

Begin 

gain:=0.0, 
loss:=0.0; 
net_benefit;=0 0; 
for k:=l to z do 

if (I[k] <= \yj) then 

gain;=gain+d[k], 
for j;=l to y do 
begin 

if (j in \yj) then 

for i;=l to x do 

if (j in mp]) then 
begin 

loss=loss+b[i], 

m[i]:=[], 

end; 

end; 

net_benefit=gain-Ioss, 

End; 
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Procedure to Find the Extent of Increase possible to Improve the Solution 


Procedure Find_extent_of_increase; 
var 

i , j , k . integer, 


Begin 

for k =1 to z do 
begin 

eoi[k] =inf, 
c22[k]:=inf; 

end; 

for k:=l to z do 
begin 

for j:=l to y do 

if not(j in wj) then 

if (c22[k]>c20,k]) then 
c22[k]:=c2[j,k]; 
eoi[k];=c22[k]-c2star[k], 
if (eoi[k]=0) then 
eoi[k]:=inf; 

end; 

eoil =inf, 
for k =1 to z do 

if (eoil>eoi[k]) then 
eoil :=eoi[k]; 
for i;=l to x do 

forj:=l toy do 
begin 

dec[ij]:=inf; 

eod[i];=inf; 



cll[i]:=inf; 


end, 

fori;=l to x do 
begin 

forj:=l toy do 

if (j in wj) then 
begin 

dec[i j] :=cl [ij]-clstar[i]; 
if (cl[ij]-clstar[i]=0) then 
dec[ij]:=inf, 

end; 

end; 

i=l; 

while(i<=x) do 
begin 

over:=true; 
eod[i]:=dec[i,l]; 
for j:=l to y do 

if (j in vvj) then 
begin 

if (dec[ij]=inf) then 
begin 

\ 

eod[i]=inf, 

over.=faIse, 

end 

else 

if (over) then 
begin 

- if (eod[i]>dec[iJ]) then 

eod[i];=dec[ij]; 


end; 



end; 


end; 

eodr=inf; 
for i:=l to x do 

if (eodl>eod[i]) then 
eodl;=eod[i]; 
if (eodl>eoil) then 
feoi:=eoil 

else 

feoi;=eodl; 


End; 

Procedure to Calculate the Best Result Obtained 

Procedure Show_Best_Result ; 
var ij,k:integer; 

Begin 

findstars; 

opt:=0.0; 

for i:=l to x do 

opt :=opt+c 1 star[i] *b [i] ; 
for k =l to z do 

opt;=opt+c2star[k] *d[k]; 

End; 



Procedure to Update the Cost Coeflicients after Finding Extent of Increase 


Procedure Updatecosts; 
var i j,k;integer, 

Begin 

find_extent_ofJncrease; 
for j;=l to y do 

if (j in wj) then 
begin 

for i:=l to x do 

cl[ij]:=cl[iJ]-feoi; 
for k;=l to x do 

c20,k]:=c20,k]+feoi; 

end; 

show_best_result. 

End; 


Main Program Begins Here 


Begin 

Initialize; 

Read_Input; 


100 : 

for k:=l to z do 
begin 

u:=u+l; 

find_stars; 

find_sets_m_andj; 

wj:=l[k]; 




on 



find_net_benefit; 
if (net_benefit>0) then 
begin 

iterations:=iterations+ 1 ; 

update_costs; 

goto 100; 

end, 

end; 


200 : 

for k:=l to z do 
begin 

u:=u+l; 

find_stars; 

find_sets_in_and_l; 

wj:=[]; 

for k2:=l to z do 
begin 

for j;=l to y do 

if (G in l[k]) and G in l[k2])) then 
Ayj:=v^+l[k2]; 

end, 

find_net_benefit, 
if (net_benefit>0) then 
begin 

iterations:=iterations+l ; 

updatecosts; 

goto 200; 

end; 

end; 

writeln('Heuristic Terminates here' ); 



write(Tinal Optimal Solution is as follows 
writeln(opt); 

writeCNo. of Iterations taken by heuristics : '); 

writeln(iterations); 

writeln; 


End. 


Main Program Ends Here 
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Problems Tested 




Problem Size : 200 x 200 x 200 
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Problems Tested 


